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Recent developments on the four dimensional (4d) lattice studies of the finite temperature electroweak phase 
transition (EWPT) are summarized. The phase diagram is given in the continuum limit. The finite temperature 
SU(2)-Higgs phase transition is of first order for Higgs-boson masses mn < 66.5 ± 1.4 GeV. Above this endpoint 
only a rapid cross-over can be seen. The full 4d result agrees completely with that of the dimensional reduction 
approximation. The Higgs-boson endpoint mass in the Standard Model (SM) would be 72.1 ± 1.4 GeV. Taking 
into account the LEP Higgs-boson mass lower bound excludes any EWPT in the SM. A one-loop calculation 
of the static potential in the SU(2)-Higgs model enables a precise comparison between lattice simulations and 
perturbative results. The most popular extension of the SM, the Minimal Supersymmetric SM (MSSM) is also 
studied on 4d lattices. 
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1. INTRODUCTION 

The visible Universe is made of matter. This 
fact is based on observations of the cosmic dif- 
fuse 7-ray background, which could be larger than 
the present limits, if boundaries between "worlds" 
and "antiworlds" had existed The observed 
baryon asymmetry of the universe was eventually 
determined at the EWPT (| . On the one hand 
this phase transition was the last instance during 
which baryon asymmetry could have been gen- 
erated around T « 100 GeV, on the other hand 
at these temperatures any B-l-L asymmetry could 
have been washed out. The possibility of baryo- 
genesis at the EWPT is a particularly attractive 
one, since the underlying physics can be -and al- 
ready largely has been- tested at collider exper- 
iments. Thus, the detailed understanding of this 
phase transition is very important. 

A succesfull baryogenesis scenario consists 
three ingredients, the Sakharov's conditions. 

1. Baryon number violating processes 

2. C and CP violation 

3. Departure from equilibrium. 

All of the three conditions has non-perturbative 
features and are studied on the lattice (e.g. at this 
conference baryon number violating sphalerons 
have been discussed by |^, spontaneous CP vi- 
olation by 1^, whereas this contribution mostly 
studies the out of equilibrium condition). 



It is rather easy to see the necessity of the 
first two conditions. Without baryon number vi- 
olation no net baryon asymmetry can be gener- 
ated. C and CP violation are needed to give a 
direction to the processes. The standard picture 
concerning the third condition is the turn-off of 
the baryon number violating rate after the phase 
transition, which means a smaller sphaleron rate 
than the Hubble rate. Inspecting the formula 
for the sphaleron rate one needs a strong enough 
phase transition, thus v/Tc>l. This ratio is of 
particular interest and both perturbative and lat- 
tice studies have the main goal to determine it. 

The first-order nature of the EWPT for light 
Higgs bosons can be shown within perturbation 
theory. However, perturbation theory breaks 
down for Higgs boson masses {ttih) larger than 
about 60 GeV due to bad infrared behavior of 
the gauge-Higgs part of the electroweak theory 
1^. Solutions of gap-equations even suggest an 
end-point scenario for the first order EWPT §. 
Numerical simulations are needed to analyze the 
nature of the transition for realistic Higgs bosons. 

One very succesfull possibility is to construct 
an effective three dimensional (3d) theory by us- 
ing dimensional reduction, which is a perturba- 
tive step. The non-perturbative study is carried 
out in this effective 3d model The end-point 
of the phase transition is determined and its uni- 



versality class is studied g] . 

Another approach is to use 4d simulations. The 
complete lattice analysis of the SM is not feasi- 
ble due to the presence of chiral fermions, how- 
ever, the infrared problems are connected only 
with the bosonic sector. These are the reasons 
why the problem is usually studied by simulat- 
ing the SU(2)-Higgs model on 4d lattices, and 
perturbative steps are used to include the U(l) 
gauge group and the fermions. Finite tempera- 
ture simulations are carried out on lattices with 
volumes ■ where Lt <C Lg are the tempo- 
ral and spatial extensions of the lattice, respec- 
tively. Systematic studies were carried out for 
TO// « 20 GeV, 35 GeV, 50 GeV and 75 GeV |]. 
The lattice spacing is basically fixed by the num- 
ber of the lattice points in the temporal direction 
(Tc = l/(Lta), where is the critical temper- 
ature in physical units); therefore huge lattices 
are needed to study the soft modes. This prob- 
lem is particularly severe for Higgs boson masses 
around the W mass, for which the phase transi- 
tion is weak and typical correlation lengths are 
much larger than the lattice spacing. In this case 
asymmetric lattice spacings are used |ll0| . 

2. END-POINT IN FOUR DIMENSIONS 

The 4-d SU(2)-Higgs model is studied on both 
symmetric and asymmetric lattices, i.e. lat- 
tices with equal or different spacings in temporal 
{at) and spatial (os) directions. The asymmetry 
of the lattice spacings is given by the asymmetry 
factor ^ = as/at- The different lattice spacings 
can be ensured by different coupling strengths in 
the action for time-like and space-like directions. 
The action reads in standard notation 

S[U,v] = 

^« E (l - ^Tr + /3t ^ (l - ^Tr C/^, 

sp ^ ^ tp ^ 

+ E|^Tr(^+^,)-f A iTr(^+^,)- 
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Figure 1. Imaginary part of first Lee- Yang zero 
as a function of A from simulations on symmetric 
lattices with Lt = 2. Filled symbols are with- 
out A-reweighting, while open symbols with A- 
reweighting from filled symbol with same shape. 
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We introduce and (3^ = (3s(3t- The 

anisotropies 7^ = A / /3s and 7k = Kt / t^s are func- 
tions of ^. We use ^ = 4.052, which corresponds 
to 7k = 4 and 7^ = 3.919. 

The determination of the end-point of the fi- 
nite temperature EWPT is done by the use 
of the Lee- Yang zeros of the partition function 
Z. Near the first order phase transition point 
the partition function reads Z = Zg + Z}, (x 
exp(— + exp(— y/h) , where the indices s(b) 
refer to the symmetric (broken) phase and / 
stands for the free-energy densities. We also 
have /{, = fs + oi[k — Kc) , , since the free- 
energy density is continuous. It follows that 
Z cx exp[-y(/s + /fc)/2] cosh[-ya(K - Kc)] , 
which shows that for complex k Z vanishes at 
Im(K) — 2TT-{n~l/2)/{Va) for integer n. Incase 
a first order phase transition is present, these Lee- 
Yang zeros move to the real axis as the volume 
goes to infinity. In case a phase transition is ab- 
sent the Lee- Yang zeros stay away from the real 
K axis. Denoting kq the lowest zero of Z^ i.e. the 
position of the zero closest to the real axis, one 
expects in the vicinity of the end-point the scal- 
ing law Im(«;o) = C[Lt,\)V-'' + Kg(Lt,A). In 
order to pin down the end-point we are looking 




Figure 2. Dependence of Rhw,c, on 1/Lj and 
extrapolation to the continuum limit. 



for a A value for which Kq vanishes. In practice 
we analytically continue Z to complex values of 
K by reweighting. Small changes in A were taken 
into account by reweighting. The dependence of 
Kq on A is shown in fig. 1. To determine the 
critical value of A i.e. the largest value, where 
Kq — 0, we have performed fits linear in A to the 
non-negative Kq values. 

In the isotropic case we have used Lt = 2. 
The Lee- Yang analysis gave Ac = 0.00116(16) for 
the end-point. Performing T = simulations 
with the same parameters this can be converted 
to rriH c = 73.3 ± 6.4GeV. In the anisotropic lat- 



tice simulation case |12 we also performed a con- 
tinuum extrapolation for Lt — 2,3,4,5 (fig. 2), 
moving along the lines of constant physics (LCP), 
and obtained 66.5 ± 1.4 GeV, which is our final 
result for the end-point in the SU(2)-Higgs model. 

Based on previous 4d simulation results one 
can determine the phase diagram of the finite 
temperature EWPT and compare it with the 3d 
analysis (fig. 3.) as it has been done in ref. 
13 1 . The phase transition lines Tc(m//), are in 
perfect agreement for mn > 25 GeV. For strong 
first order phase transition close to the Coleman- 
Weinberg limit the 3d approach seems to be less 
accurate. The error bars on the endpoints are on 
the few percent level, thus uncertainty of the di- 




Figure 3. A comparison of the phase diagrams 
obtained from direct 4d simulations (squares) 
and from dimensionally reduced 3d simulations 
(shaded region). 



mensional reduction around the end-point is also 
in this range. This indicates that the analogous 
perturbative inclusion of the fermions results also 
in few percent error on mn- 

One can determine what is the endpoint value 
in the full SM. As it was shown previously the 
perturbative integration of the heavy modes is 
correct within our error bars. Therefore we use 
perturbation theory to transform the SU(2)-Higgs 
model endpoint value to the full SM. We obtain 
72.1 ± 1.4 GeV, where the dominant error comes 
from the measured error of RHW.cont.- The error 
on gj^ is eliminated by calculating the relationship 
between the coupling definitions used in pertur- 
bation theory (MS) and lattice simulations (from 
static potential) ||l^,|lj| . The calculation of this 
relationship and a comparison of the perturbative 
and lattice results on the EWPT will be shortly 
discussed in the next section. 

The full SM result needs some explanation. 
Based on vacuum stability the measured top mass 
{nitop ~ 175 GeV) results in a lower bound for the 
Higgs boson mass (approx. 130 GeV). This value 
is higher than the previously mentioned 72.1±1.4 
GeV. For the pure SU(2)-IIiggs model the end- 
point Higgs mass is 66.5 ± 1.4 GeV. The inclusion 



of the fermions, especially the top increases the 
endpoint slightly. For a hypothetical top quark 
mass less than approximately 150 GeV the lower 
bound is less than « 70 GeV, thus it is below 
the endpoint and it gives a reliable theory. In- 
creasing the top quark mass the lower bound gets 
larger than the endpoint. This means that inde- 
pendently of the direct experimental bounds on 
the Higgs boson mass no EWPT exists in the SM. 

3. RELATIONSHIP BETWEEN GAUGE 
COUPLINGS 

Despite the fact that the perturbative and lat- 
tice approaches are systematic and well-defined, 
it is not easy to compare their predictions. The 
reason is that in lattice simulations the gauge cou- 
pling constant is determined from the static po- 
tential, whereas in perturbation theory the MS 
scheme is used. One can calculate the static po- 
tential on the one-loop level in the SU(2)-Higgs 
model 13,1^]. As expected the numerical dif- 
ference between the two conventions is not that 
large, it is within a few percent, for details see 
p3| , p^ . With this connection we could perform 
a precise comparison between the predictions of 
perturbative and lattice approaches (fig. 4). 

In jlj] the existing lattice data was reanalyzed 
and a continuum limit extrapolation was per- 
formed whenever it was possible. The only quan- 
tity which is measured so precisely that the def- 
inition of the gauge coupling constant is essen- 
tial is the ratio of the critical temperature to the 
Higgs boson mass. As it has been observed al- 
ready for Mh ~ 35 GeV the perturbative value of 
Tc is larger than in lattice simulations. This sort 
of discrepancy disappears for larger Higgs boson 
masses. A plausible reason for this fact is the con- 
vergence of the high temperature expansion used 
in the perturbative approach. 

The most dramatic differences appear clearly 
as we get closer to the end point. The pertur- 
bative approach gives non-vanishing jump of the 
order parameter, non-vanishing latent heat and 
interface tension, while the lattice results suggest 
rapid decrease of these quantities and no phase 
transition beyond the end-point. 
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Figure 4. Pulls plotted against the Higgs mass 
for Tc, w/T, latent heat and interface tension. 

4. PHASE TRANSITION IN THE MSSM 

As it was demonstrated in the previous sec- 
tions the SM is not suitable for baryogenesis, 
not even for a first order EWPT. Several ex- 
tended models were studied in order to obtain 
a stronger first order phase transition and a reli- 
able baryon asymmetry. The most popular model 
is the MSSM, which perturbatively shows a much 
stronger phase transition than the SM (even 
an intermediate colour breaking phase transition 
is possible in these scenarios) . Lattice studies 
in a 3d reduced model (with one Higgs doublet) 
basically confirmed the perturbative results ||l^ . 

We performed a 4d lattice study with the 
bosonic sector of the MSSM Q. The lattice ac- 
tion is too long to be presented here, thus only 
the fields involved are listed. Both of the Higgs 
doublets, the stop, sbottom scalars and SU(2), 
SU(3) gauge fields were included. It is of particu- 
lar importance to keep both of the Higgs doublets, 
since according to the standard scenario the gen- 
erated baryon asymmetry is directly proportional 
to the change of the ratio of their expectation val- 
ues ub oc {v'^)A(3{Tc). Here the length squared 
of the Higgs field (w^ = -|-z;|) is integrated over 
the bubble wall. The ratio of the expectation val- 
ues of the two Higgs fields is tan/? = Vi/v2, and 
the difference between the /3 values are taken in 
the "symmetric" and in the "broken" phases. 

We had simulations at Lt = 2, 3, 4, 5 and moved 
along the line of constant physics. Our simulation 
point corresponds to tan/3(r = 0) « 6, and the 
mass of the lightest Higgs bosons is approx. 35 
GeV (in the bosonic theory). 

Two values of as were taken (the physical 
and a smaller one). The physical as resulted in 



v/Tc « 1.5, whereas the smaher value of gave 
a stronger phase transition v/Tc ~ 2. Perturba- 
tion theory predicts just the opposite behaviour 
(stop-gluon setting sun graphs are proportional to 
the strong coupling and they are responsible for 
the strengthening of the phase transition). The 
reason can be the difference between the renor- 
malization effects in the stop sector. 

We measured the (3 parameter in both phases at 
the phase transition. One obtains tan^ j3{sym) — 
38.13(15), iim^l3{hrok) = 36.04(15), which gives 
A/3 = 0.0045(7). This result is far below the per- 
turbative prediction A(3{pert.) — 0.017. 



5. CONCLUSIONS 

The endpoint of hot EWPT with the technique 
of Lee- Yang zeros from simulations in 4d SU(2)- 
Higgs model was determined. The phase tran- 
sition is first order for Higgs masses less than 
66. 5 ± 1.4 GeV, while for larger Higgs masses only 
a rapid cross-over is expected. The phase diagram 
of the model was given. 

It was shown non-perturbatively that for the 
bosonic sector of the SM the dimensional reduc- 
tion procedure works within a few percent. This 
indicates that the analogous perturbative inclu- 
sion of the fcrmionic sector results also in few 
percent error. In the full SM we get 72.1 ± 1.4 
GeV for the end-point, which is below the lower 
experimental bound. This fact is a clear sign for 
physics beyond the SM. 

Based on a one-loop calculation on the static 
potential of the SU(2)-Higgs model a direct com- 
parison between the perturbative and lattice re- 
sults was performed. 

The MSSM is more promising for a succesfuU 
baryogenesis. Some 4d results were shown, indi- 
cating a strong first order phase transition. 
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